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Abstract 

Randall-Sundrum models based on SU (2) l x SU (2)^ with custodial symmetry are com- 
pelling frameworks for building alternative models of electroweak symmetry breaking. 
A particular feature of these models is the likely presence of light Kaluza-Klein fermions 
related to the right-handed top quark. These can be as light as a few hundred GeV and 
still compatible with EW precision constraints. In this article, we study the detectability 
of four-TV final states at the LHC, which arise from the pair-production and tW decay of 
light Kaluza-Klein bottom quarks as well as light Kaluza-Klein quarks carrying electric 
charge 5/3. 



1 Introduction 



In the last eight years, extra-dimensional models have been suggested to solve the gauge 
hierarchy problem. Among them, the Randall-Sundrum (RS) model [1] is the most appealing, 
where the hierarchy between the electroweak (EW) and the Planck scales arises from a warped 
higher dimensional spacetime. Variants of the original set-up have matured over the years. 
Eventually, all Standard Model (SM) fields except the Higgs (to solve the hierarchy problem, it 
is sufficient that just the Higgs -or alternative dynamics responsible for electroweak symmetry 
breaking- be localized at the TeV brane) have been promoted to bulk fields rather than brane 
fields. It has been shown that EW precision constraints are much ameliorated if the EW 
gauge symmetry in the 5-dimensional bulk is enlarged to SU{2) L x SU(2) R x U(l)x [2]. The 
AdS/CFT correspondence suggests that this model is dual to a strongly coupled CFT Higgs 
sector [3]. Also, the SU(2) L x SU(2) R gauge symmetry in the RS bulk implies the presence of 
a global custodial isospin symmetry of the CFT Higgs sector, thus protecting EW observables 
from excessive new contributions [2] . This gauge structure in warped space has also been used 
to construct Higgsless models of EW symmetry breaking [4] . 

In this framework, Kaluza-Klein (KK) excitations of gauge bosons of mass Mkk ~ 3 TeV 
are allowed (even lower, ~ 500 GeV, for Higgsless models). Generically, the mass spectrum 
of fermionic KK modes depends on the boundary conditions (BC) imposed on the TeV and 
Planck branes, localized at the boundaries of the slice of AdS^. These are commonly modelled 
by either Neumann (+) or Dirichlet (— ) BC in orbifold compactifications. Early studies in 
Randall-Sundrum background considered Neumann boundary conditions only, in which case 
all KK masses have to be larger than the gauge boson KK mass of 3 TeV and will be difficult 
to produce at the LHC. In recent models, different boundary conditions have been imposed, 
leading to richer possibilities for model building and to the potential presence of light KK 
fermions, observable at colliders. These fields have for instance Dirichlet BC on the Planck 
brane and Neumann BC on the TeV brane. They are denoted ( — h) KK fermions, and do not 
have zero modes. 

The heaviness of the top quark is explained by the localization of the wave function of 
the top quark zero mode near the TeV brane, guaranteeing a large Yukawa coupling with the 
Higgs. In the initial model of Ref. [2], the Right-Handed (RH) top quark is included in a 
doublet of the SU(2) R symmetry. Its b R partner does not have a zero mode but its first KK 
excitation is expected to be light. This mode mixes with the SM bottom quark and induces 
large corrections to Z — > bb, so its mass has to be at least ~ 1.5 TeV [2]. However, it has 
been recently pointed out in Ref. [5] that the custodial symmetry, together with a discrete 
L <-> R symmetry and alternative SU(2) R assignments for the top and bottom quarks, can 
protect the Zbb coupling and allow light masses for the KK b R as low as a few hundreds of 
GeV. This mode is accompanied by other light degenerate KK quarks (carrying electric charge 
Q = 2/3,-1/3,5/3), that have been named "the custodians" [6]. They are likely to be the 
lightest KK states in these models and could be produced at the LHC. 

The interesting phenomenology of light ( — h) KK fermions was pointed out in Ref. [7] and 
its extended version Ref. [8], where it was emphasized that light KK fermions are expected 
as a consequence of the heaviness of the top quark. More precisely, KK partners of the 
RH top should be light. Even though the study of [7,8] deals with the embedding of these 
SU(2) L x SU(2) R xU(l) Randall-Sundrum models into a GUT0, the statement remains pretty 

1 In this case, an extra Z3 symmetry can be imposed to protect proton stability, leading to the stability of 
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general: ( — h) KK partners belonging to the tn multiplet are expected to be light and can 
be probed at colliders. The possibility of pair-producing these KK fermions, leading to multi 
W events, was discussed. Most of these signatures were specific to the GUT model except 
for the production of the KK (that we denote which is common to a large class of 
models. In this paper, we focus on the pair production and tW decay of bn, which leads to a 
very unique AW + bb final state. This is a promising signature for probing RS models with 
custodial symmetry^!. Moreover, while we will focus on the b~R in this paper, it should be noted 
that our study applies as well to pair production of other custodians such as the KK b^ and 
the exotic KK quark with electric charge Q = 5/3, that both decay into tW. 

The pair production of exotic quarks, singlet under SU(2) L , was considered in the past 
in Ref. [10]. The single and pair production of heavy T quarks in Little Higgs models was 
studied in Ref. [11]. However, these quarks have charge +2/3 instead of -1/3 or 5/3, so the 
decay is into Wb rather than Wt and there are only 2 W's in the event. We are not aware 
of previous studies investigating events with more than 2 Ws in the final statfl While A-W 
production brings the prospect of a spectacularly rich multi-jet plus multi-lepton signature, 
distinguishing the signal from backgrounds will still be challenging. Our goal in this paper is 
to investigate the feasibility of doing so and thus identifying bn at the LHC. 

A straightforward trigger criterion for these events is that of a single, isolated lepton with 
missing Et, i.e., the standard leptonic W data stream. Other W's in the event would be 
reconstructed using dijet pairs. Multi-lepton triggers can also be used. The fermions bn are 
particularly accessible, producing b jets which could be tagged with displaced tracks. If one 
bu enters the leptonic W data stream, its antiparticle partner is potentially amenable to full 
reconstruction. 

In the present work, we neglect hadronization effects. For the range of masses considered, 
the width of bn is large enough that the decay takes place before hadronization. 

For low bft masses (< 300 GeV), there are interconnection effects. This situation is common 
to nearly all investigations of heavy particles with GeV-scale widths such as WW, ZZ, and 
it where hadronic decays overlap. We assume this effect is small; this is justified for the 
(9(500) GeV masses we consider in detail, and we neglect it for now (it is interesting to note 
that on the opposite extreme, the KK fermions of [8] are very long-lived and lead to CHAMP- 
like signatures; for these, we are well into the realm of Heavy Quark Effective Theory, so that 
to first order we can again neglect light quark effects in the decay). 

In Section [21 we describe the parameters of the model and the decay channels of b^. 
Readers who are more keen on the experimental aspects than on the model-building details 
can skip this section. Section [3] discusses the possible signatures associated with b^ pair 
production and presents the main SM backgrounds mimicking the AW signature. Section [4] 
gathers the results of our simulation and outlines a promising strategy for distinguishing the 
signal from the SM backgrounds using early (10 fb _ ) LHC data. 

a stable light KK right-handed neutrino. 

2 Collider signatures of this class of models are only starting to be investigated, see for instance [9]. 

3 Multiple weak gauge boson production was theoretically investigated in the eighties [12, 13] and early 
nineties [14-16] in the context of the Superconducting Super Collider (SSC), as these processes were (and are) 
of great interest to test whether the EW symmetry breaking is due to new strong interactions. If this is the 
case, multiple (> 2) longitudinal W and Z production would be enhanced at SSC energies [12, 14]. This is 
nevertheless still a small effect compared to the rates for the standard production of multiple transverse W 
and Z and it was concluded that it would not be observable [13]. More recent work has suggested that these 
studies could be performed with 0(1000) fb _1 of data that could be delivered by an upgraded LHC [17]. 
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2 Model parameters 



Our analysis describes the pair production and decay of b R . However, as we said in the intro- 
duction, there are typically other "custodian" quarks which lead to the same AW signature, 
potentially significantly enlarging our signal. The number of custodians is model-dependent. 
Following [5], several papers have appeared recently on the phenomenology associated with 
these light KK fermions [6,18-20]. These studies consider different embeddings for the SM 
top and bottom quarks, leading to various possibilities for the number and masses of the 
custodians. In the description below, we assume for simplicity that there is only one light 
KK quark, b R , the SU(2) R doublet partner of the SM RH top quark. This is not the realistic 
situation but it does not matter for our analysis. At the end of this section, we comment how 
to generalize our discussion to the more realistic models mentioned above. 

KK fermions are 4-component spinors which have both a LH and a RH chirality. One 
chirality has (— +) BC while the other has (-1 — ). The RH chirality of the KK b R under 
consideration turns out to be localized near the TeV brane and the LH one near the Planck 
brane. Therefore, the couplings of the LH chirality with modes localized near the TeV brane, 
such as the Higgs, the top quark and KK excitations, are suppressed. In contrast, the direct 
interactions of KK fermions to zero-mode gauge bosons are vector-like. This is because zero- 
mode gauge bosons have a flat profile (unlike KK modes) and couple identically to both 
chiralities of KK fermions. The couplings of the KK b R to SM gauge bosons are the same as 
the couplings of the SM b R quark, except that they involve both chiralities identically. The 
vertices &r&rG m , b R b R A^ and b R b R Z^ are respectively (? S 7 M , —e^^/S and e tan ^7^/3- bR has 
four decay channels: b R — > t R W, b R — > tj^W, b R — > b L H and b R — > b^Z. The first two lead 
to the multi W signature we are interested in, while the third can also lead to it if the Higgs 
decays dominantly into WW. We now present these decays in details. 




2.1 b R — > t R W decay 

b R , as a singlet under SU(2)l, is not expected to couple to W. The effective coupling of b R to 
W is due to W R — Wl mixing resulting from EW symmetry breaking (see Fig. [3]). W R is the 
charged gauge boson associated with SU{2) R and does not have a zero mode. We consider 
the effect of the mixing of the SM W with the first KK excitation of W R , which has mass 
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Figure 2: The four decay channels of &r. 



Mkk ~ 3 TeV. The induced coupling is then g~ b t w _ 7^ where: 



9R I 

9l R ,t R ,w- = -7K\/ k7rr c x P R x M Wr -w l x 2\ 



b R ,t R 



(1) 



M Wr -w l = —\jlkixr c 




(2) 



Ti R t ~ 1 is the form factor reflecting the overlap between the wave functions of 6r, W# and 
£r. M.y/ R -w L is the mixing factor due to the EW breaking vev of the Higgs. knr c = ln(10 15 ) is 
the exponent of the warp factor needed to generate the weak/Planck hierarchy. This formula 
assumes the Higgs is localized on the TeV brane. If it is delocalized in the extra dimension 
as in scenarios of gauge-Higgs unification where it is identified with the 5th component of a 
gauge boson, we can replace the factor \/2/c7rr c with the form factor accounting for the wave 
function overlap. In this case, using the profile given in [21] (see also Appendix B of [8]) 
\J1kixr c ~ 8.31 is replaced by 5.748. Pr is the RH projector that expresses the fact that only 
one chirality of the KK has a non-suppressed coupling to Wr. gR is the 4D gauge coupling 
of SU(2)r. The LR discrete symmetry requires = gR at 5D level. Since the 5D gi is fixed 
by the matching with the 4D gi (assuming that brane kinetic terms are small), we can replace 
9r by 9-, the SM SU (2) coupling, in the above formulae. As is clear from Fig. 0], this decay is 
negligible compared to the other three channels. 



2.2 bu —> tiW , bR — > biZ and bR — > b^H decays 

The Yukawa coupling between the (i^, hi) and (tR, bR) multiplets generates three more decay 
modes for bR-. bR — > b^H , bR — > t L W~ and b r — > b^Z. 



\ 1 
\ / 




Figure 3: Diagram leading to W R — mixing. 
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The coupling between bi, H and 6r is y t f(c tR ) /\/2, where y t = 1 and f(c) ~ J2/ (1 — 2c) 
(for c > —1/2) comes from the wavefunction of the RH top quark given in appendix A 
of [8]. The c-parameter corresponds to the 5D fermionic bulk mass in Planck units. For the 
right-handed top quark, we take c = —1/2, leading to f(c tR ) = 1. 

The second decay comes from the interaction bR.t[ J .H ± where H is the would-be Gold- 
stone boson which becomes a longitudinal W. In unitary gauge, there is no H , and this has 
to be written in the form of a 7 M interaction. However, such an interaction will couple a LH 
(RH) fermion with a LH (RH) fermion. As we said earlier, the RH chirality of the KK b R 
under consideration turns out to be localized near the TeV brane and the LH one near the 
Planck brane. Therefore, the coupling of the LH chirality with the top quark is suppressed 
(since the top is localized near the TeV brane). One would be tempted to conclude that there 
is no coupling between bn ti and W, but there is a subtlety: the LH chirality of the physical 
(mass eigenstate) bn is actually a linear combination of a LH KK bf> and a zero mode bi- 
There is some mixing effect (comparable to the one discussed in section 9.3 of [8]) due to the 
top Yukawa coupling which generates a mass term between b^ and 6r. After diagonalization 
of the mass matrix, we can see that the new mass eigenstate 5r has an admixture of the zero 
mode through which it couples to W and Z. Specifically, in unitary gauge, 

[physical b] LH = cos 9 b R + sin6> 6 L (3) 

where 5r denotes the LH chirality of KK bn, and sin^ « m t f(c tR )/m~ bR comes from the mass 
term between b L and Sr. Via the bi component, the physical b couples to W and t L (with SM 
coupling g/y/2). When we project onto the longitudinal W, we get a factor of E/mw from 
the polarization vector, where E ~ m l R smce we are considering decays of 5r. The coupling 
to Wi ong is y t f(c tR ) as expected using the goldstone equivalence theorem (coupling to charged 
Higgs). To summarize, in unitary gauge, the physical b R has the following interactions: 
b R {g Z b L b L Z + {g/y/2) t L W)m t f(c tR )/m~ bR where g Z b L = (g/ cos6' H /)(l/2 - sin 2 8 w /3). 

2.3 Decay widths 

The partial width for the decay 6r — » t^W - is 



dwtr 2 



m bR + ( m w 2 ~ 2m t 2 )m- b 2 + m t 4 - 2m w 4 + m t 2 m w 2 
?t R w- = ^ 2 A [vm,m w ,mi J (4) 

where 

/ f m w \ 2 , w \ i , ( m t- m w) 2 ( m t+ m w) /c\ 

9wtr = 9r kirr c — — and \{m t , m w , mi = 1 H 2 = . (5) 

\M KK J R ml R mf R 

We use the same formula for the decay 6r — > tiW~ except that the coupling constant is 
replaced by 

ft \ ^ 9 /«X 

9wu = f(c tR )—-x- (6) 

The decay width for bn — > b^Z is 

9zbi 2 (m- bR 2 - m z 2 ) (m^ 4 + m z 2 m^ 2 - 2m z A ) 
bLZ ~ 32vr ml m z 2 [ ) 

OR 
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where 

9m = f{c tR )— x - 9 [3 - 2 sin 2 9 W ] 
m- h 6 cos &w 



and finally 



bLH 327rm? 



2 

(9) 



"R 



As we said, in our analysis, we use f(c tR ) = 1, g R = g, Mkk = 3 TeV, and kiir c = ln(10 15 ), 
so that the only free parameter is rn~ bR . 

In our formulae, we have assumed that (t#, 6r) form an SU(2)r doublet even though we 
know that this is not the realistic situation. The solution to the Zbb problem indeed requires 
that tR is either a singlet under SU (2) R or belongs to a (1,3) + (3, 1) under SU (2) L x SU(2) R 
[5]. In the first case, there is obviously no b R . However, this does not mean that there is no 
associated light KK quark. Indeed, in gauge-Higgs unification models (see Ref. [21]) which are 
presumably the best motivated models for EW symmetry breaking in RS, the SM t R belongs to 
a larger multiplet which necessarily leads to light custodian partners. In minimal composite 
Higgs models, t R belongs to a 5=((2,2),(1,1)) or a 10=((2,2),(1,3) + (3,1)) of SO(5). Its 
(2, 2) partners contain the SM Ql as well as an extra doublet of custodians of electric charge 
(2/3,5/3). The Q=2/3 quark will not lead to the AW signature, but the Q=5/3 quark (that 
we name q) will, and with a branching ratio essentially equal to 1. If t R belongs to a 10, there 
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will be as well two bottom- like custodians b^ and b^, and two alike q custodians all leading 
to the AW signature. All these particles have the same mass, so they are produced with the 
same cross section. The L «-> R discrete symmetry also guarantees the same couplings for 
bft and bi and for q^ and q~i- However, the fact that they belong to a triplet rather than a 
doublet will change the couplings presented below by a Clebsh-Gordan factor of order 1. In 
any case, for our analysis, the total widths are not important: the 6r, bi and q decay almost 
as soon as they are produced. 

In our simulation, we have chosen to illustrate the signal in the case where tji belongs 
to (1,3) + (3,1). The number of AW events at the end will thus be 2(1 + B 2 )/B 2 times the 
number of events obtained from just one b^ decaying to tW (as given in Fig. [6]), where B is 
the relevant branching ratio. 

To summarize this discussion, whatever model is being considered, there will be at least one 
light KK quark, related to the SM top quark, that will lead to our AW signature. Additional 
light KK quarks will also contribute. Note also that the value of the Higgs mass depends in 
principle on the choice of embedding for the top and bottom quarks. In Ref. [6,18], ran and 
m~ bR are strongly related. Models with gauge-Higgs unification typically have a light Higgs, 
while for the present analysis we allow the Higgs mass to be free and as heavy as 300 GeV. 

3 Preferred decays and event selection 

In this study, we focus on those decay channels which show promise for being detected in the 
pp collisions of the LHC First, we rely on a single isolated lepton arising from one of the W's 
to provide a clean and efficient trigger. We then reconstruct other W's in the event using 
pairs of jets. W bosons mostly decay into a pair of quarks (BR = 67.96%), while they decay 
into a charged and neutral lepton pair about 10.68% of the time for each lepton generation. 

The actual rate of multi-jy events arising from our signal depends upon the bn mass as 
well as the Higgs mass: for instance, in the case of mi < 500 GeV, the branching ratio into 
bZ increases dramatically from the ~ 20% characteristic of higher b^ masses, leading to decay 
signatures such as 2 b jets + missing Ej-, 2 b jets + 2 leptons + missing E?, or 2 b jets + 4 
leptons. We do not pursue these signatures here. 

For mi R > 500 GeV, however, the Higgs mass plays an important role. If ran ~ 115 
GeV, it decays mainly into b quarks, and the b^ — > bH channel leads to 6 b jets. Since the 
branching ratio for bH is typically ~ 50%, this takes half of the decays away. On the other 
hand, if the Higgs decays mainly into WW, then even the 6r decay into bH produces the 4 
W + 2 b signature. In this case, most of the produced b r lead to this signature. Note that 
we are assuming that b^ is the lightest of the KK fermions, as this is the natural situation 
in minimal models. However, it could happen, such as in the GUT models of Ref. [7,8], that 
there are lighter KK states such as the KK RH neutrino. The Higgs would then mainly decay 
into this KK RH neutrino, the H — > WW branching ratio would practically vanish, and we 
would lose the possibility of increasing the WW yield from the Higgses. 

If the Higgs mass is in the region where the dominant decay is actually to WW*, the 
visible signal rate will also be reduced, since W* — > jj will not peak very much around the W 
mass. Of course, if W* —>l~v and W — > jj, we still get our trigger and our W combination. 
For ran > 300 GeV, H — > WW — > Ivjj becomes significant, so ran = 300 GeV is a suitable 
choice for a first look at multi-lV signatures. The branching fraction of H — > WW is about 
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70% in the m# region of 200 to 350 GeV, where the ti channel opens up. 

It should also be noted that in a typical detector, it is not trivial to measure the charge 
of a jet's precursor; moreover, the jet energy is reconstructed with a finite resolution. For 
instance, for dijet pairs with pt > 350 GeV, the ATLAS Technical Design Report [22] claims 
a mass resolution of 6.9 GeV, which is not much less than the difference between the W and 
Z masses. The two mass peaks will overlap significantly, and indeed may be seen as a single 
bump. The proposed detection method will therefore have difficulty differentiating genuinely 
WWWW events from, for instance, WWZZ arising from a decaying to bH, followed by 
H — > ZZ. For the purpose of discovering bn, it is not clear that there is any advantage 
in discriminating between W's and Z's here. In the present signal simulation, we do not 
consider these possibilities, but in principle, we should simulate these other decays as they 
can potentially increase our signal. 

3.1 Backgrounds 

The Standard Model background for AW production (see Ref. [15,16]) is suppressed by high 
powers of electroweak couplings and we neglect it. The dominant background is actually 
fake and comes from the misinterpretation of jets as coming from W's. The two important 
backgrounds we consider arise from ti and tiH production, ti leads to 2 W's and 2 6's with 
four extra jets misinterpreted as coming from hadronic decays of W's. At large Higgs mass, 
where the Higgs decays primarily into WW, tiH exactly leads to AW + 2-6-jets and becomes 
a serious background for heavy b^ masses: if m~ bR ~ 1 TeV, the two production cross sections 
are comparable [23]. Distinguishing the signal in this case would be a significant challenge. 

Diboson EW production can also contribute by misidentification and event overlaps. How- 
ever, the LHC production cross sections for these processes are smaller with respect to ti 
production. In addition, they do not produce a pair of high p t b-jet pairs. Therefore, contri- 
butions from such processes are not currently estimated. A list of diboson production cross 
sections taken from Ref. [24] are listed below along with that of it as taken from the ATLAS 
Technical Design Report and references therein [22,25]: 



The cross section for gg — ► tiWW is 4 orders of magnitude below that of tt and this 
background is neglected. Similarly, the AW + Ab background from (ti) (ti) production can be 
neglected as it is more than 4 orders of magnitude smaller than the ti cross section: the phase 
space volume is scaled down by a l/(47r) 4 factor and there is a g\ suppression factor in the 
cross section. In the end, this is a 2 s /(A-K) 2 suppression compared to the ti cross section. This 
is consistent with the cross section obtained with Comphep: 1.3 x 10 -2 pb just from the gg 
contribution. 

Finally, there are of course Standard Model EW processes which give W + 6-jet events, 
but these are of higher order, and in addition they have no b jets so we ignore them. There 
are also QCD processes which will give a hard photon + 6 jets, with the photon converting 
to a high-pr lepton. Such events should not have much missing E t and so should also be 
suppressed. 



tt 

W+W 

wz 

ZZ 



833 pb (NLO + NLL resummation) 



86.7 pb 
32.4 pb 
12.9 pb 
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4 Signal and background simulation 



The Lagrangian of the model is implemented into CalcHEP (v2.4.3 [26]), a tree level Monte 
Carlo event generator that can deal with multi-particle final states. We generate tiWW events 
from bR pair production in CalcHEP, which are further processed in PYTHIA 6.4.01 [27]. 




, I I I , I I I , I I , , I , , , I , , , I , , , I , I , I I I I I , •» I 

200 400 600 800 1000 1200 1400 1600 1800 2000 
M(B R ) (GeV) 



Figure 5: Pair production cross section for bn at the LHC. Only masses above the t + W 
threshold are considered. 

Figure [5] shows the bft pair production cross section at the LHC. We use CTEQ6L parton 
distribution functions and consider only the contribution from light quarks, u and d. The 
dominant contribution comes from the gluon-initiated hard scattering up to about a 5r mass 
of 800 GeV. The QCD scale was set to the mass of Sr. The cross sections include the 
contributions from s-channel EW exchange diagrams, though the combined contribution from 
these channels is negligibly small. For example, for 6r mass of 500 GeV, the contribution 
from EW channels is about three orders of magnitude smaller than those involving the QCD 
coupling. 

Figure [6] shows the 4W + 26 final states expected in lOfb^ 1 of LHC data as a function 
of the bfi mass, for two Higgs masses and considering only the production cross section and 
branching fractions. For a 6r mass of 500 GeV, the yields are about 5000 events from tW and 
900 from bH for a Higgs mass of 300 GeV. The yield is about 2800 events from only the tW 
channel for a Higgs mass of 115 GeV, with no contribution from the bH channel. 

As discussed in Section 12. 3[ there are additional exotic quarks with exactly the same 
properties as 6r which therefore enhance the total number of events: for rn~ bR = 500 GeV and 
ran = 300 GeV, including b^L and qn t L, this multiplying factor amounts to 11.2. In the case 
of the models discussed in Ref. [6,18-20], where some have only one light Q = 5/3 quark q 
instead of &r, this factor is 1/ B 2 = 4.6. 
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* tW (m H =300 GeV) 
tW(m H =115GeV) 
*• bH (m H =300 GeV) 
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200 400 600 800 1000 1200 1400 1600 1800 2000 
M(b R ) (GeV) 



Figure 6: Number of 4W + bb events expected from the tW and bH decays of bft for 10 fb 1 
of LHC data. 

4.1 Trigger strategy and event reconstruction 

In order to investigate how such events could be examined using an LHC detector such as 
ATLAS, we first restrict ourselves to simulating 6r pair production and their decays through 
the tW channel. For the present model analysis, we choose mi R = 500 GeV and m# = 300 
GeV. We apply event filters and acceptance criteria based on nominal ATLAS parameters as 
found in the Technical Design Report [22]. The following "trigger", applied to the generated 
events, is based on the lepton criteria for selecting W — > Iv events: at least one electron or 
muon with pt > 25 GeV must be found within the pseudorapidity range \r)\ < 2.4, where 
rj = — lntan(0/2) and 9 is the angle relative to the pp collision axis; then, the "missing 
Et" , calculated by adding all the neutrino momenta in the event and taking the component 
transverse to the collision axis, must exceed 20 GeV. 

We mimic cone-based hadronic jets as they might be observed in a detector: stable charged 
and neutral particles within \rj\ < 4.9 (the range of the ATLAS hadronic calorimeter), exclud- 
ing neutrinos, are first ranked in order. Jets are seeded starting with the highest pt tracks, 
with pt > 1 GeV; softer tracks are added to the nearest existing jet, as long as they are 
within AR < 0.4 of the jet centroid, where AR = \/ Acf) 2 + Arj 2 . The number of jets with 
Pt > 20 GeV is shown in Figured The signal is peaked around 8 jets. Additional jets can be 
produced in quark hadronization and underlying parton activity, while jets can be lost due to 
falling below the energy threshold or outside the geometric acceptance. 

The background sample is dominated by it events generated using TopRex (version 4.11) 
[28] and PYTHIA (version 6.403), with CTEQ6L parton distribution functions. The small 
itH contribution to the background has been modelled with PYTHIA. As expected, the 
background has fewer high-p^ jets than the signal, peaking around 5 jets. 

Figures M and M compare signal and background distributions for transverse momenta of 
generated particles as well as event-level observables such as the "trigger" (highest-pr) lepton 
Pt, missing Et, and scalar sum of the Et's of all the jets in the event. 
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Figure 7: Number of jets with > 20 GeV for signal and background. Both distributions 
are normalized to unit area. 
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Figure 8: Generator-level signal and background distributions after the "trigger" condi- 
tions, normalized to unit area. Top left: t quarks. Top right: b quarks. Bottom left: W 
bosons. Bottom right: W daughter quarks. 
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Figure 9: Generator- level signal and background "trigger" distributions after initial trigger 
conditions, normalized to unit area. Top left: Pt of stiffest electron or muon. Top right: 
missing Et- Bottom left: scalar sum of Et in the event. Bottom right: number of W's 
decaying hadronically in the event. 
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The scalar sum distribution in Figure [9] (bottom left) suggests that the background 
can be significantly reduced by requiring the sum to exceed 800 GeV. We also restrict the 
analysis to events with between 6 and 9 jets. No b jet tagging is performed at this point. 

After applying the 800 GeV cut to the scalar Et sum, we construct hadronic W candidates 
by adding the 4-momenta of two jets, each with > 20 GeV and assuming that each 
individual jet has zero mass. Since the W's are typically produced with p T > 150 GeV, as 
shown in Figure [SJ the total pt of the dijet combination is further required to exceed that 
value. The mass spectrum for the signal and background events, representing 10 fb _1 of LHC 
data, is shown in Figure [10J We have scaled up the signal to include the and q relatives 
of the b R , as discussed in Section 12.31 (this scaling assumes charge-symmetric lepton and jet 
identification). A prominent peak can be seen around 80 GeV with width 5 GeV. The width is 
due to, among other effects, the hadronization of the daughter quarks, with the resulting jets 
sometimes overlapping with other activity in the event. The low mass of the peak, relative to 
the generated W mass, is expected given, for instance, the finite cone size of the jets. A W 
peak is also evident in the tt-dominated background distribution, as expected. Figure [10] also 
shows the reduction in the signal's combinatorial background if perfect b tagging is achieved, 
though experimentally, b tagging is rarely very efficient or pure. 




Figure 10: Left: dijet mass distribution for all jets for signal, background, and their sum, for 
10 fb _1 integrated luminosity. Right: comparison of dijet mass distributions before and after 
eliminating all b jets. The black histogram here is the same as the blue in the figure to the 
left. 

In order to suppress the most common (tt) Standard Model background, we eliminate 
the single hadronic W decay in the following manner: we start with the highest px jet and 
search among the lower p? jets for a combination whose mass falls between 70 and 90 GeV. 
If no combination is found, the search is continued using the next highest px jet. If a pair 
is found, those two jets are removed for consideration in forming the subsequent dijet mass 
combinations, which are shown in Figure [TTJ The observable peak is by now dominated by 
the signal, though the background peak has not been entirely eliminated. These background 
events might appear, paradoxically, to contain 3 W's, one decaying leptonically and two 
hadronically, but a more mundane explanation is that the trigger lepton and missing energy 
actually arise from semileptonic b decay occurring among the far more numerous tt events. 
These events also appear in Figure [9] (bottom right) as background events with 2 hadronic W 
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decays; a similar fraction of signal events contain 4 hadronic W decays. While the background 
peak may not be very large compared to the signal, it can still be further suppressed by 
increasing the trigger thresholds. 




Figure 11: Dijet mass distribution and signal, background, and their sum, after eliminating 
the first hadronic W candidate. The distribution represents 10 ftT 1 integrated luminosity. 

Relating the 5r production cross section to the size of the peak (and therefore the discovery 
potential of this signature) is non-trivial: first, we would have to subtract or suppress the tt 
component, model the dijet mass resolution, and calculate the reconstruction efficiency in light 
of the multiple W decays in the event. Such details are best reserved for further investigation, 
when a more detailed and specific simulation of the detector and experimental environment, 
including the effects of detector material and multiple pp interactions in a beam crossing, can 
be employed. Furthermore, the relationship between 6r production and yield is complicated 
by the m^-dependent interplay of different decay modes and how they contribute W's and Z's 
to the signal. However, this initial analysis, indicative of the kind which could be attempted 
at an LHC experiment, suggests that it will be possible to distinguish the signal from the 
major physics backgrounds. 

4.2 Further analysis directions 

It has been suggested in Ref. [18] that if all four Vy's decay leptonically, the resulting "golden" 
signature of 4 leptons + bb + missing Et would have very little SM background. As shown 
in Figure [12], however, the rate for this signature would be suppressed by approximately two 
orders of magnitude relative to a signature requiring only one leptonic W decay. At the same 
time, because of the multiple neutrinos, it would be difficult to identify the leptons as W 
daughters. As a result, while this signature could be indicative of new physics, it would be 
difficult to distinguish what kind of new physics is being observed. 

An important piece of corroborative evidence for these models is the observation of the 
custodian quark q carrying electric charge Q — 5/3. Pair production of this quark would lead 
to exactly the same 4 W + bb final state as the 6r. A possible method to distinguish the 
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Figure 12: Number of 4W events (coming from b^ only) in the case ra# = 300 GeV, with the 
assumption of 10 fb _1 integrated luminosity. Leptons are either e or \i. 

q events is to select events with two high-p^ leptons with the same sign (presumably from 
the two W's from a single q), and then fully reconstruct the other q in the event through 
its tW — > bWW — > bjjjj decay. This method is driven by the typical situation that the 
charge of an observed lepton can be measured reliably, while that for a jet cannot. The 
fully reconstructed q would yield a narrow peak in the b + 4 jet combined invariant mass 
distribution, while the corresponding distribution arising from bf> pair production and decay 
would be much broader. The main disadvantage of this method is the typically low efficiency 
of a full reconstruction, which requires all decay products to be observed and well measured. 
We therefore expect that confirming the existence of the q custodian would require rather 
more data than the 10 fb~ studied here (note, however, that the LHC's design luminosity 
would yield 100 fb _1 per year). 

Finally, let us note that one could also consider the single production of 5r together with 
the standard b quark via s-channel Z or W^-gluon fusion, or together with the top quark 
via s-channel W. This would be especially relevant for large b r masses and would deserve a 
separate study. 

5 Conclusion and future prospects 

We have studied new signals in pair production of heavy Q = —1/3 and Q = 5/3 quarks at the 
LHC. They are produced through standard QCD interactions with a cross section ~ O(10) 
pb for masses of several hundreds of GeV. Heavy quarks such as b^ are well-motivated in 
Randall-Sundrum models with custodial symmetry. They are generally Kaluza-Klein part- 
ners of the Standard Model Right-Handed top quark. Their decay channels were described in 
detail in this paper. In the present work, we focussed on the 4-W events which we believe are 
quite specific to this class of models, and also experimentally promising. We have considered 
the process gg, qq — > b R b R — > W~t W + t — > W~W + b W + W~b where at least one W boson 
decays leptonically and the other ones hadronically. A simulation of this signal and its main 
background was performed, and an analysis strategy outlined which distinguishes the signal 
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from the sizable Standard Model backgrounds. The peak we obtain in the dijet mass distri- 
bution suggests that it is possible to reach a signal significance beyond the 5a level. Further 
study with more detailed simulation is required to map the discovery potential for this signal 
at an LHC experiment such as ATLAS, or at the ILC, and to connect the observable signal 
to the production cross section. 
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